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ABSTRACT 


A  method  is  described  for  calculating  the  activation  energies  of  bimolecular 
reactions  of  such  multivalent  gaseous  compounds  as  O^,  N^,  NO,  N^O,  NO^, 
CO,  and  CO^.  It  is  also  applicable  to  the  monatomic  reactants  O,  N,  and 
H  and  to  the  OH  radical.  The  method  uses  relationships  between  dissocia¬ 
tion  energy,  bond  length,  vibrational  wavenumber,  and  bond  order.  It  avoids 
the  use  of  adjustable  parameters.  The  agreement  between  calculated  and 
experimental  activation  energies  is  very  good  when  spin  conservation  in  the 
transition  state  is  considered  in  the  computation  procedure.  A  convenient 
rule  has  also  been  formulated  that  predicts  the  activation  energy  will  be- 
negligible  in  bimolecular  reactions  between  radicals  when  the  transition- 
state  formula  corresponds  to  that  of  a  compound  having  stxble  bonds.  No 
exceptions  to  this  rule  have  been  found.  The  method  is  applicable  to  com¬ 
puting  rates  needed  for  nonequilibrium  calculations  of  propellant  performance 
and  of  reentry  behavior. 
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I.  INTRODUCTION 


The  reaction  rates  of  compounds  containing  oxygen,  nitrogen,  or  car¬ 
bon  are  important  in  determining  the  nonequilibrium  effects  on  propellant 
performance  of  the  fuels  and  oxidizers  composed  of  these  elements.  In 
reentry  through  the  earth's  atmosphere,  particularly  if  ablation  of  carbon¬ 
aceous  materials  is  involved,  the  rates  of  reactions  of  such  compounds  must 
affect  reentry  nonequilibrium  phenomena.  These  rates  can  also  be  of  con¬ 
siderable  significance  in  affecting  propellant  ignition,  jet  flames,  super¬ 
sonic  wind  tunnels  and  upper-astmosphere  composition.  They  can  also  in¬ 
fluence  areas  of  peripheral  aerospace  interest  such  as  air  pollution  and  the 
internal  combustion  engine. 

Although  the  activation  energies  and  rate  constants  for  some  of  the 
reactions  of  these  compounds  have  been  measured,  a  basis  for  predicting 
and  understanding  their  kinetic  behavior  has  not  been  established.  For  ex¬ 
ample,  the  activation  energies  of  the  exothermic  reactions  NO^  +  O  =  NO  + 
and  NO^  +  H  =  NO  +  OH  are  zero,  but  the  activation  energies  of  the  two 
exothermic  reactions  +  NO  =  N-,  +  NO.,  and  O ^  +  CO  =  O  +  CC^  are  each 

50  kcal/mole.  This  difference  of  50  kcal/mole  in  activation  energy  provides 

35 

the  huge  factor  of  more  than  10  toward  reducing  the  rate  constant  of  the 
latter  two  reactions.  Since  the  preexponential  factor  of  the  Arrhenius  rate 
constant  equation  is  closely  related  physically  to  the  readily  calculated  colli¬ 
sion  frequency,  the  preexponential  factor  can  usually  be  predicted  within  a 
factor  of  less  than  5  for  similar  reactions.  Consequently,  the  most  impor 
tant  problem  in  the  prediction  of  rate  constants  for  these  multivalent  reac¬ 
tions  is  that  of  computing  the  activation  energies.  The  computational  method 
developed  in  this  study  provides  an  explanation  for  these  highly  significant 
differences  in  activation  energies.  The  method  is  also  to  be  applied  in  the 
computation  of  rate  constants  of  the  important  O,  N,  and  C  compounds  in  those 
reactions  for  which  no  experimental  kinetic  data  exist. 
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Transition- state  computations  ’  of  activation  energies  and  rate  con¬ 
stants  for  hydrogen-atom  transfers  in  the  gas  phase  exhibit  good  agreement 
with  measured  activation  energies  and  rate  constants  when  several  trial 

assumptions  involving  bond  energy,  bond  order,  and  bond  length  interrela- 

3 

tionships  are  used.  Included  are  the  assumptions:  (1)  that  the  Pauling 
"rule"  for  the  relationship  between  bond  length  and  bond  order  applies  in 
the  transition  state;  (2)  that  the  Lennard- Jones  6-12  energy  and  distance 

4 

parameters  for  the  noble  gases  represent  the  zero  bond  order  energies 
and  internuclear  distances  for  elements  in  the  corresponding  row?  of  the 
periodic  tables;  (3)  that  an  approximately  linear  relationship  exists  between 
the  logarithm  of  bond  energy  and  the  logarithm  of  bond  order*;  (4)  that  the 
triplet  repulsion  energy  arising  from  uncoupled  spins  in  the  transition  state 
is  represented  by  a  Sato'’  anti-Morse  function  similar  to  that  for  triplet 
Instate  of  H£.  Instead  of  having  to  use  these  assumptions  for  computing 
kinetic  parameters,  a  quantitative  quantum-mechanical  treatment  would 
clearly  be  preferable, but  a  tractable  quantum  treatment  does  not  exist  for 
complex  multielectron  reacting  systems.  It  is,  nevertheless,  possible 
(although  open  to  serious  question  )  that  the  changes  in  potential  energy  of 
dissociating  bonds  or  forming  bonds  may  be  described,  as  a  function  of  bond 
length  or  nominal  bond  order,  within  some  arbitrary  error  such  as  2  kcal 
per  bond  by  postulates  derived  empirically  from  energy  data  for  a  series 
of  bonds  covering  a  range  of  bond  orders.  ^ 

The  aforementioned  activation  energy  computing  method  does  not 
utilize  adjustable  parameters.  It  does,  however,  use  a  tentative  model  of 
the  transfer  reaction 

AH  +  B  =  A-..  H--B  =  A  +  HB  (1) 

ni  n2 

whereby  the  sum  of  the  bond  orders  rij  and  n^  in  the  transition  state  is  con¬ 
stant  at  unity,  corresponding  to  the  initial  univalent  bond  order  of  the  reac¬ 
tant  AH  and  of  the  final  bond  order  of  the  product  HB.  For  this  model,  the 
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activation  energy  of  the  transition  state  is  calculated  for  one  line,' 
rather  than  a  surface  as  in  the  London- Eyring-Polanyi  method^*  of  attempt¬ 
ing  to  compute  activation  energy. 

12 

Because  of  the  relative  success  *  of  this  bond-energy  bond-order 
(BEBO)  model  in  predicting  rate  constants  for  H-atom  transfers,  it  was 
decided  to  investigate  its  applicability  to  reactions  in  which  the  transferred 
univalent  atom  was  larger  than  the  hydrogen  atom  and,  eventually,  to  reac¬ 
tions  in  which  multivalent  atoms  were  transferred.  Thus,  propellant  and 
reentry  reactions  involving  com  pounds  of  O,  N,  F,  Cl,  C,  etc.,  could  be 

7 

covered.  A  reduced- variable  treatment  was  introduced  so  that  larger, 
univalent  atoms  such  as  the  halogens  were  scaled  relative  to  hydrogen.  The 
computed  activation  energies  and  rate  constants  were  found  to  agree  very 
well  with  the  measured  values  of  those  relatively  few  halogen  transfer  reac¬ 
tions  for  which  the  experimental  kinetic  data  were  available.  The  number 
of  such  reactions  was  too  few  to  produce  much  confidence  in  the  modified 
BEBO  method  when  not  limited  to  H-atom  transfers,  particularly  since  the 
measured  activation  energies  were  all  quite  low  (3  kcal/mole  or  less). 

There  is,  however,  a  relatively  large  quantity  of  kinetic  data  available 
for  atom  transfer  reactions  involving  oxygen  or  nitrogen  atoms,  and  the 
measured  activation  energy  barriers  for  these  transfers  include  rather  high 
values  as  well  as  low  values.  Modifications  of  the  method*  were  needed 
because  of  the  multivalent  character  of  the  dissociating  bonds  or  forming 
bonds  involving  nitrogen  or  oxygen  in  the  transition  state.  The  multivalency 
clearly  voided  the  use  of  the  simple  BEBO  reaction  model  in  which  it  was 
postulated  that  n^  equals  1  -  n  for  univalent  bonds  in  the  transition  state. 

This  assumption  is  equivalent  to  stating  that  the  bond  order  of  H*  •  •  B 
is  equal  to  the  decrease  in  bond  order  n  j  of  the  dissociating  bond  A*  ♦  •  H 
from  its  initial  bond  order.  An  extension  of  this  trial  assumption  to 
multivalent  transfers  can  be  made  by  choosing  a  transition-state  model 
in  which  the  bond  order  n.,  continues  to  be  linearly  related  to  the  decrease 
in  the  bond  order  of  the  dissociating  bond  A*  •  •  X,  where  X  is  the  trans¬ 
ferred  atom.  The  computation  procedure  associated  with  this  model  is 
described  in  the  next  section.  Comparison  of  experimentally  determined 
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activation  energies  with  those  computed  from  this  model  with  spin  con¬ 
servation  applied  to  the  transition  state  showed  very  good  agreement, 
when  it  is  considered  that  the  formula  of  the  transition  state  in  some 
cases  corresponded  with  the  formula  of  a  stable  compound  (e.  g. ,  . 

Activation  energies  were  computed  for  other  bond  dissociation -bond  forma¬ 
tion  models  such  as  those  corresponding  to  nonconservation  of  spin, or  bond 
energies  rather  than  dissociation  energies.  The  computation  procedures 
for  these  other,  relatively  unsatisfactory  models  are  described  in  the  latter 
portions  of  the  next  section. 
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II.  COMPUTATION  PROCEDURES 


A.  MULTIVALENT,  DISSOCIATION  ENERGY 


The  aforementioned  transition- state  model  for  multivalent  transfers 
in  which  the  bond  order  of  the  formed  bond  X*  •  •  B  was  linearly  dependent 
on  the  decrease  in  bond  order  of  the  dissociating  bond  was  satisfied  by  the 
relationship 


n2  =  <nXB/nAX*  <nAX"nl) 


(2) 


where  nAX  and  n^g  are,  respectively,  the  initial  bond  order  of  the  reactant 
AX  and  the  final  bond  order  of  the  product  XB,  X  is  the  transferred  atom 
(e.g. ,  nitrogen  or  oxygen),  and  n^  is  the  bond  order  of  A»**X  in  the  transi¬ 
tion  state.  Equation  (2)  meets  the  requirement  that  when  n^  is  zero  (at  the 
initiation  of  the  reaction)  n^  equals  n^.and  that  n^  equals  n^g  when  n^  is 
zero  at  the  conclusion  of  the  bimolecular  transfer  reaction.  When  X  is  a 
univalent  atom  such  as  hydrogen,  Eq.  (2)  reduces  to  the  BEBO  trial  assump¬ 
tion  that  nj  +  ng  equals  1. 

In  order  to  compute  activation  energies  for  transfers  involving  a 
multivalent  bond  by  a  modified  BEBO  method,  it  was  necessary  to  extend  the 
bond-energy  bond-order  relationship  to  orders  above  1.  Since  a  linear 
relationship  between  log  Dg  (where  D#  is  the  gas-phase  dissociation  potential 
energy  of  the  bond)  and  log  n  was  used  for  bond  order  not  exceeding  1,  a 
relatively  simple  modification  of  the  BEBO  method  was  made  by  also  using 
a  linear  log  D^/log  n  relationship  between  n^  (or  n^g)  and  n  equal  to  1. 
Then, the  slope  of  the  line  for  the  multivalent  region  is  given  by 


«AX  1  t  log  (D^/Dj)  3/log  t,AX 


where  is  the  dissociation  potential  energy  of  the  bond  A-X  of  bond  order 
1.  D  includes  the  sero-point  energy,  as  shown  in  Table  1,  since  it  is  based 

*  g 

on  the  minimum  of  the  potential  energy  curve.  For  diatomic  species,  the 
dissociation  energy  is  equal  to  the  bond  energy  B^,  but,  as  shown  in  Table  I, 


TABLE  I.  Molecular  and  Bond  Data3, 


Data 

Molecules 

Nz 

NO 

°2 

CO 

no2 

n2° 

N(NO) 

cT1 

O 

Dissociation  energy, 

D  ,  kcal/ mole 

228.6 

152.8 

120.3 

258.9 

73.7 

41.0 

116.2 

128.4 

t  -8 

Bond  length,  rg,  10  cm 

1.098 

1.151 

1.207 

1.128 

1.  197 

1. 186 

1.126 

1. 159 

Vibrational  wavenumber, 

2358 

1904 

1580 

2170 

1330 

1410 

1685 

1720 

UJ  .  cm'1 
e 

Bond  energy,  B^,  kcal/mole 

228.  6 

152.8 

120.3 

258.9 

113.2 

110 

158 

193.5 

Bond  order 

3.0 

2.5 

2.0 

3.0 

1.75 

1.75 

2.5 

2.0 

Ground  state 

L 

2* 

3Z 

4 

2Z 

4 

4 

Single-bond  values  used 

Dissociation  energy, 

39.7 

49.5 

35.0 

85.7 

32.2 

18.  1 

29.1 

55.6 

D,  ,  kcal/mole 

Bond  length,  rj, cm 

Vibrational  wavenumber, 
cm'1 

1.48 

1.44 

1.47 

1.43 

1.44 

1.44 

1.48 

1.43 

1370 

1220 

1120 

1180 

1080 

1090 

1105 

1180 

Bond  ener;  ,  Bj,  kcal/mole 

39.7 

49.  5 

35.0 

85.7 

49.5 

49.5 

39.7 

84.0 

a  r he  JANAF  Thermochemical  Tables  (with  supplements  to  December  1966)  were  the  primary 
source  of  data.  Energies  were  computed  on  the  basis  of  the  potential  energy  curve  and,  conse¬ 
quent!  ,  ‘ncluded  the  zero-point  energy.  In  polyatomic  molecules,  De  for  the  bond  broken  or 
formed  is  not  equal  to  the  bond  energy  since  the  other  bonds  in  the  molecule  change  in  energy 
v  hen  one  bond  is  broken  or  formed.  For  NgO,  which  can  dissociate  to  form  either  Njr  or  NO,  the 
total  gas-phase  dissociation  energy  of  268  was  partitioned  between  the  ON  and  NN  bonds  on  the 
basis  of  their  estimated  bond  orders. 
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this  equality  does  not  hold  for  polyatomic  species  since  the  bond  energies  of 

the  other  bonds  in  the  molecule  generally  change  when  a  bond  is  broken  or 
3 

formed.  The  values  of  D  in  Table  I  were  obtained  in  a  straightforward 

e  8 

manner  from  thermochemical  data,  but  the  values  bf  Dj  for  the  triatomic 
molecules  were  estimated  because  they  have  not  been  measured.  The  estima¬ 
tions  were  based  on  the  simple  assumption  that  the  relationship  of  Dj  to 
was  in  the  same  proportion  as  D  to  B 


D./B.  =  D  /B 
11  e  e 


(4) 


This  expression  also  implies  that  the  slope  *s  unchanged  if  B^Bj  is  used 
in  Eq.  (3)  instead  of  D  /D,  . 

With  these  modifications,  the  potential  energy  V  of  transition- state  for¬ 
mation  in  multivalent  bimolecular  transfer  reactions  was  calculated  in  a  man- 

1  2 

ner  similar  to  that  used  for  univalent  transfers.  * 


V  =  <De.AX-Dl,AX"l 


where  Vr  is  the  triplet  repulsion  as  calculated  by  the  reduced  variable  treat¬ 
ment  of  Eq.  (10)  in  Rsf.  7.  When  a  bond  order,  or  n^,  is  less  than  1,  the 

slope  n  in  Eq.  (5)  is  replaced  by  the  slope  p  calculated  for  the  bond-order 

1  2 

region  below  1.  In  a  manner  similar  to  that  of  previous  investigations,  ’ 
the  computer  program  modified  for  multivalent  bonds  determines  the  value 
of  n 2  that  corresponds  to  the  peak  in  the  potential  energy  V  of  the  linear 
transition- state  model  A*  •  •  X  •  •  -B.  The  computed  activation  energy  for  this 
model  is  equal  to  the  potential  energy  peak  height  plus  a  small  zero-point 
energy  correction.  ^ 

B.  SPIN  CONSERVATION  IN  THE  TRANSITION  STATE 

The  dissociation  energies  Dg  shown  in  Table  I  correspond,  as  is  cus¬ 
tomary,  to  dissociation  to  ground-state  products.  In  thermal  dissociation, 
the  ground-state  atomic  oxygen  dissociated  from  CO^(  I)  and  woul  ’ 
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3  1 

be  triplet  0(  P),  but  spin  would  be  conserved  if  a  singlet  0(  D)  atom  were 

dissociated  from  the  parent  molecule.  For  these  bimolecular  transfers,  it 

is  not  known  whether  dissociation  of  CO^  or  starts  by  splitting  off  a 

singlet  O  (^D)  atom  and  is  followed  by  spin  conversion  deexcitation  of 

3 

the  singlet  to  the  triplet  0(  P) .  It  is  conceivable  that  prior  to  the  dissocia¬ 
tion  the  molecule  transforms  to  an  excited  state  via  an  internal  spin  change 

3 

and  that  0(  P)  is  then  formed  directly  by  dissociating  from  an  excited  mole¬ 
cule,  or  the  oxygen  atom  dissociated  from  ground  state  or  CO^  could 
cross  over  to  O(^P).  From  a  kinetics  point  of  view,  it  appears  that  the  most 
direct  mode  of  beginning  to  break  the  oxygen  bond  in  bimolecular  transfers 
involving  N^O  or  CO^  could  well  be  by  the  formation  of  singlet  0(XD).  From 
that  point  of  view,  it  is  clearly  desirable  to  consider  a  transition -state  model, 
when  N^O  or  CO_  are  involved  in  the  bimolecular  transfer  of  an  oxygen  atom, 

CL  j 

in  which  the  dissociating  oxygen  (i.  e. ,  X  in  A •  •  •  X*  •  •  B)  is  in  the  C>(  D)  state. 

This  can  be  done  quite  readily  in  computations  of  the  activation  energy  by 

3  1 

adding  the  difference  in  energy  between  0(  P)  and  0(  D)  (45.4  kcal/mole)  to 

D  when  the  described  computation  procedure  is  used  for  multivalent  transfers, 
e 

In  contrast  with  the  dissociation  spin-conservation  requirements  for  N_0 

3  c 

and  CC^,  it  is  possible  for  ground  state  0(  P)  to  dissociate  directly  with 
spin  conservation  from  the  other  molecules  of  Table  I:  NO,  CO,  NO^,  and 

CL.  Consequently,  a  transition-state  model  for  oxygen  transfers  involving 

^  .  1 
only  the  latter  molecules  need  not  consider  the  possibility  that  0(  D)  is  the 

state  of  the  transferred  c  tom. 

C.  BOND  FORMATION  IN  THE  TRANSITION  STATE 
WITHOUT  CORRESPONDING  DISSOCIATION 

12  7 

When  the  modified  BEBO  method  was  applied  ’  ’  to  univalent  bonds 
such  as  those  in  or  H^,  the  model  of  the  transition  state  logically  involved 
the  assumption  that  dissociation  of  the  bond  of  reactant  AX  must  occur  in 
proportion  to  the  formation  of  the  product  bond  XB.  This  proportionality  is 
expected  since  the  XB  bond  is  formed  with  the  participation  of  an  electron 
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that  had  been  used  in  the  AX  bond.  In  multivalent  reactions,  the  situation  can 

l 

be  more  complex.  For  example,  in  the  case  of  NC^  reactant,  the  valence  elec¬ 
trons  are  not  in  a  completely  bonded  form,  as  indicated  by  the  electron  spin. 

It  is  possible,  therefore,  that  a  model  for  the  oxygen-transfer  reaction  of  NO^ 
with  H  could  consist  of  the  formation  of  a  species  such  as  H"  *ONO,  where  the 
bond  order  of  H  •  •  -O  increases  without  a  similar  decrease  in  the  ON  bond.  In 
this  study,  the  criterion  for  selecting  those  transition  states  for  which  this 
model  was  considered  possibly  applicable  was  the  existence  of  a  stable  compound 
that  corresponded  to  the  transition- state  formula.  In  the  cited  case,  nitrous 
acid  would  be  the  stable  compound,  with  a  dissociation  potential  energy  of  78 

g 

kcal/mole  for  the  HO  bond.  Other  compounds  with  formulas  similar  to  those 
of  relevant  transition- state  models  are  N^O,  NC^,  and  O^.  Computations  of 
activation  energy  based  on  this  model  can  be  made  by  subtracting  a  term  from 
the  right-hand  member  of  Eq.  (5)  that  corresponds  to  the  potential  energy  of 
the  bond  formed  by  one  reactant  with  the  other  reactant  (e.  g. ,  H  and  ONO). 

D.  BOND- ENERGY  DISSOCIATION 

The  dissociation  models  previously  used*  for  polyatomic  molecules  allow 
the  dissociation  of  an  atom  from  a  reactant  to  proceed  slowly  enough  so  that 
the  other  bonds  in  the  reactant  reach  their  ground  state  before  the  A*  •  •  X  bond 
(where  A  is  not  monatomic)  is  completely  broken.  Consequently,  the  dissocia¬ 
tion  energy  rather  than  the  bond  energy  was  used  for  computations  based  on 
this  model.  It  is  conceivable,  however,  that  for  some  polyatomic  molecules 
the  bond  A*  •  *X  may  break  so  rapidly  in  the  transition  state  that  the  other 
bonds  in  A  remain  unchanged  in  the  transition  state.  When  computations  were 
made  for  oxygen  or  nitrogen  transfers  based  on  this  model,  the  calculated 
activation  energies  did  not  usually  agree  well  with  the  measured  values.  For 
this  model,  computations  were  made  by  using  the  bond  energy  data  of  Table  I 
rather  than  the  dissociation  energies. 
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E.  SIMULTANEOUS  POLYVALENT  BOND  FORMATION 

AND  DISSOCIATION 

It  is  characteristic  of  the  transition- state  models  previously  des¬ 
cribed  that  bond  formation  proceeds  first  via  formation  of  the  single  bond, 
which  then  becomes  a  polyvalent  bond  as  the  reaction  progresses.  It  is 
clearly  possible  that  the  formation  of  a  polyvalent  bond  in  the  transition 
state  does  not  proceed  through  any  recognizable  single-bond  state  such  as 
those  for  which  data  are  presented  in  Table  I, but  that  the  molecular  orbital 
changes  involved  in  the  formation  of  a  polyvalent  bond  proceed  directly 
from  the  zero  bond  to  the  final  bond  without  the  intermediate  step  of  passing 
through  a  bond  characteristic  of  single  bonds  in  stable  molecules.  Similarly, 
the  dissociation  of  the  polyvalent  bond  A*  •  «X  in  the  transition  state  could 
proceed  directly  to  the  zero-bond  state  without  following  the  potential-energy 

bond-order  path  suggested  by  the  dissociation  energy  D.  of  a  single  bond. 

1  1 

By  using  the  D  and  r  of  Table  I  and  the  Johnston-Parr  trial  assumption 
e  e  ^ 

that  the  Lennard- Jones  6-12  energy  and  distance  parameters  for  neon  rep¬ 
resent  the  zero  bond-order  energies  and  internuclear  distances  of  N  and  O 
bonds,  values  of  and  were  calculated  for  use  in  Eq.  (5)  over  the 
entire  bond-order  range  rather  than  limiting  the  use  of  an<*  ffXB  to  k°nc* 
orders  above  1.  The  activation  energies  computed  for  this  model  gave  poor 
agreement  with  experiment  for  most  of  the  reactions  studied,  but  not  for 
reactions  of  an  atom  with  a  diatomic  species.  In  the  latter  cases,  the  transi¬ 
tion  state  is  triatomic,  and  the  only  bonds  in  the  transition  state  are  the 
breaking  and  forming  bonds. 
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in.  RESULTS  AND  DISCUSSION 


Computed  activation  energies  are  compared  in  Table  II  with  experi¬ 
mental  activation  energies;  none  of  the  reactions  in  Table  n  are  of  the  atomic 
plusdiatomic  type.  The  computed  Eq's  of  particular  interest  are  in  the  spin 
conservation  column.  This  column  summarizes  the  computed  Eq's  obtained 
by  combining  the  postulate  of  spin  conservation  for  dissociation  of  the  multi¬ 
valent  bond,  described  in  the  section  on  Computation  Procedures,  with  the 
linear  bond-order  relationships  for  multivalent  bonds  as  embodied  in  Eqs. 

(2)  through  (5).  The  agreement  between  the  Eq's  computed  by  this  method 

and  the  experimental  Ert's  is  usually  quite  good.  Among  this  first  group  of 

u  1  . 
reactions,  spin  conservation  involves  the  formation  or  reaction  of  an  0(  D) 

atom  if  CO^  or  N^O  is  a  reactant  or  product,  and  the  agreement  between 

computed  and  measured  Eq's  is  much  superior  when  spin  conservation  is 

assumed.  The  multivalent  dissociation  energy  column  summarizes  the  Eq's 

computed  with  the  spin  conservation  ignored.  The  latter  dissociation  model 

provides  satisfactory  agreement  only  for  those  reactions  in  which  the  dissocia 

tion  and  formation  of  the  bonds  can  be  accomplished  via  ground- state  species 

as  in  the  NO  +  NO  =  N  +  NC^  reaction. 

The  no  triplet  column  presents  Eq's  computed  by  the  procedure  used 
for  the  spin  conservation  column,  except  that  triplet  repulsion  was 
arbitrarily  set  at  0.  From  the  results  listed  ir.  this  column  it  can  be  seen 
that,  for  almost  all  of  these  reactions,  the  triplet  repulsion  made  only  a 
small  contribution  to  the  activation  energy,  relative  to  the  contributions  of 
the  bonding  terms  in  Eq.  (5).  The  single-bend  energy  column  lists  Eq's 
computed  by  the  procedures  used  for  the  spin  conservation  column,  except 
that  the  single-bond  energy  Bj  was  used  for  Eqs.  (3)  and  (5),  rather  than  the 
corresponding  dissociation  energy  Dj.  It  is  evidently  incorrect,  as  antici¬ 
pated,  to  use  Bj  rather  than  D^  since  the  computed  Eq’s  do  not  generally 
exhibit  satisfactory  agreement  with  experiment  when  Bj  is  thus  employed. 
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TABLE  II.  Comparison  of  Computed  Activation  Energies  With  Observed 

Activation  Energies. 


Reaction 

AX+B  =  A+BX 

Observed 

Eo 

Computed  Activation  Energies3 

Spin 

conser¬ 

vation 

Multi. 

dissoc. 

energy 

No 

trip¬ 

let 

Single  - 
bond 
energy 

Mult.  - 

bond 

energy 

Poly. 

simul¬ 

taneous 

n2o+no  =  n2+no2 

a 

o 

48 

2 

46 

31 

13 

1 

n2o+o  =  n2+o2 

27d 

29 

0 

26 

35 

- 

N-jOtH^  =  N2+OH 

1 5e 

15 

0 

14 

1 

0 

o2+co  =  o+co2 

50l 

49 

37 

48 

1 

73 

5 

HO+CO  r  H+C02 

0* 

1 

0 

0 

- 

- 

no2+co  =  nchco2 

32h 

28 

2 

26 

34 

0 

ONN+O  =  ON+  NO 

271,  j 

29 

24 

25 

48 

40 

0 

N02+N  =  ON+NO 

od 

- 

0 

0 

0 

1 

0 

^Energies  expressed  in  kcal/mole. 

bF.  Kaufman  and  J.  R.  Kelso,  J.  Chem.  Phys.  23,  602  (1955). 

CE.  S.  Flshburne  and  R.  Eds*.  J.  Chem.  Phys.  41,  1297  (1964). 

dD.  Gutman,  R.  L.  Belford,  A.  J.  Hay,  and  R.  Panirov,  J.  Phys.  Chsm,  70,  1793  (1966). 

Dixon-Lewis,  M.  M,  Sutton,  and  A,  Williams,  Tenth  International  Symposium  on 
Combustion,  Combustion  Institute,  Pittsburgh,  Pa.""{l965J, "p.  "49S.  — 1 — c - 

fK.  G.  P.  Sulzmann,  B.  F.  Myers,  and  E.  R.  BarUe,  J.  Chem.  Phys.  42,  3969  (1965). 

•g.  Dixon-Lewis.  W.  E.  Wilson,  and  A.  A.  Weetenberg,  J.  Chem.  Phye.  44,  2877  (1966). 

hH.  S.  Johnston,  W.  A.  Bonner,  and  D.  J.  Wilson,  J.  Chem.  Phye.  26,  1002  (1957). 

lK.  Schofield.  Planet.  Space  Sci.  15,  643  (1967), 

■*R.  M.  Fristrom  and  A.  A.  Weetenberg,  Flame  Structure.  McGraw-ftill  Book  Co.. 

New  York,  (1965),  p.  355,  37|,  372.  - ~ 
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For  the  multivalent  bond-energy  column,  the  method  described  in  the 
Computation  Procedures  section  under  "Bond-Energy  Dissociation"  was 
used.  Examination  of  the  computed  Eq's  in  this  column  of  Table  II  shows 
that  they  do  not  generally  agree  with  the  experimental  Eq's, although  there 
are  some  cases  of  agreement.  For  the  polyvalent  simultaneous  column, 
the  method  used  was  described  in  the  Computation  Procedures  section  under 
"Simultaneous  Polyvalent  Bond  Formation  and  Dissociation."  For  those 
reactions  having  a  transition  state  with  more  than  two  bonds,  this  model  for 
the  bonding  changes  during  a  reaction  and  does  not  give  agreement  with  ex¬ 
periment. 

The  first  six  reactions  of  Table  III  are  the  only  ones  in  the  tables  for 
which  the  structural  formula  of  the  transition  state  is  similar  to  that  of  a 

Q 

stable  species  (e.  g. ,  HONO,  ^O,  NC^,  NO^,  O^,  and  HC^).  A  stable 
species  can  be  defined,  in  this  context,  as  one  in  which  no  bond  dissociation 
energy  is  less  than  25  kcal/mole.  For  this  specific  class  of  reactions,  Eq 
can  be  computed  as  outlined  under  "Bond  Formation  in  the  Transition  State 
without  Corresponding  Dissociation"  of  the  Computation  Procedures  section. 
For  this  group  of  reactions  between  two  radicals,  EQ  in  the  exothermic  di¬ 
rection  is  almost  zero,  apparently  because  sufficient  energy  is  initially  pro¬ 
vided  by  the  preliminary  partial  formation  of  the  new  bond  in  the  transition 
state  to  supply  the  subsequent  dissociation  energy  requirement  of  the  old 
bond  with  virtually  no  potential  energy  of  activation.  It  is  expected  that 
the  absence  of  activation  energy  barriers  will  be  generally  characteristic 
of  bimolecular  transfer  reactions  for  which  the  transition-state  structural 
formula  corresponds  to  that  of  a  stable  compound,  if  the  bond  formed  be¬ 
tween  the  two  radical  reactants  has  a  binding  energy  larger  than  the  activation 
energy  predicted  by  the  appropriate  model  of  Table  III.  Of  the  six  reactions 
in  this  group  in  Table  III,  four  are  the  monatomic-plus-diatomic  polyvalent 
reaction  type  for  which  the  simultaneous  bond-formation  model  may  be  justi¬ 
fied.  The  potential  value  of  the  simultaneous  model  is  further  indicated  by 
the  results  in  Table  III  for  the  last  reaction,  NO  +  O  =  N  +  O^,  which  does 
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TABLE  III.  Reactions  That  Have  a  Transition- State  Formula 
Similar  to  a  Stable- Molecule  Formula 


Reaction 

AX  +  B  =  A  +  BX 

Transition- 
state  formula 

Observed 

E0 

Computed  Eq 

_  b 

Energy 

of  new  bond 
(kcal/mole) 

Multi. 

dissoc. 

energy 

Poly. 

simul¬ 

taneous 

ON  OH  =  ON  +  OH 

hn°2 

oc 

1 

- 

78 

ON  O+O  -  0N  +  02 

NO} 

od 

3 

- 

49 

ON+N  =  0+N2 

n2o 

od 

25 

0 

61 

ON  +  O  =  O+NO 

no2 

~oe 

25 

1 

72 

o2+o  =  o+o2 

°3 

~oe 

51 

4 

26 

HO+O  =  H+O, 
c 

ho2 

of 

25 

5 

63 

Oz+N  =  O+NO 

.. 

o‘ 

27 

8 

— 

aIn  all  these  reactions,  spin  conservation  in  the  transition  state  can  occur 
without  involving  excited  atoms,  such  as  0(*D)  or  N(*D),  Therefore,  the  E. 
computed  by  the  multivalent  dissociation-energy  model  is  equal  in  these  cases 
to  the  Eq  computed  for  the  spin-conservation  model. 

**  Bond-energy  data  based  on  JANAF  tables.  In  each  of  these  six  reactions,  this 
bond  energy  is  easily  sufficient  to  supply  all  the  activation  energy  that  would  have 
been  needed,  based  on  Eq  computed  by  the  appropriate  model.  Consequently, 
observed  Eq  should  be  negligible  (i.e, ,  0  or  1  kcal/mole). 

e  t„  F,  Phillip*  and  H.  I.  Schiff.  J.  Chem.  Phys.  22  ****(1*62). 

*Y.  F.  Phillips  and  H.  I.  Schiff.  J.  Chem.  Phys,  36,  1509(1962). 

rJ.  T.  Herron  and  F.  S.  Klein.  J.  Chem.  Phys.  40,  2751(1964). 

F.  S.  Klein  and  J.  T.  Herron,  J.  Chem.  Phys.  TT,  1265  (1964). 

fR.  M.  Fristrom  end  A.  A.  Westenberg.  Flame  Structure.  McGraw-Hill  Rook 
Co..  Inc.,  New  York  (1965),  p.  155.  171.  ill, 
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not  have  a  transition- state  formula  that  corresponds  to  a  stable  compound. 
The  Eq  computed  by  the  simultaneous  model  for  this  reaction  agrees  very 
well  with  experiment,  whereas  the  Eq's  computed  by  models  involving  typical 
single-bond  energies  for  N-O  and  0-0  are  much  too  high. 

The  agreement  between  computed  and  experimental  Eq's,  particularly 
for  spin-conservation  dissociation  in  the  transition  state,  suggests  that  it 
may  now  be  possible  to  understand  important  and  previously  unrecognized 
characteristics  of  bimolecular  chemical  reactions  to  an  extent  that  was  not 
hitherto  possible.  To  summarize,  the  trial  postulates  found  to  be  of  value  in 
these  computations  are:  (1)  spin  conservation  for  the  transferred  atom  is 
maintained  in  the  transition  state;  (2)  the  activation  energy  (exothermic 
direction)  is  generally  almost  zero  when  the  structural  formula  of  the  transi¬ 
tion  state  corresponds  to  that  of  a  reasonably  stable  molecule;  (3)  the  depen¬ 
dence  of  potential  energy  on  bond  order  in  the  transition  state  is  represented 
within  about  2  kcal/mole  by  Eq.  (5),  as  supplemented  by  Eqs.  (2)  through  (4); 
and  (4)  in  reactions  of  atoms  with  diatomics  involving  polyvalent  bonds,  the 
direct,  simultaneous  formation  and  dissociation  of  polyvalent  bonds  occur. 
These  trial  postulates  provide  a  hitherto  nonexistent  basis  for  calculating 
activation  energies  and  estimating  rate  constants  of  multivalent  transfer 
reactions  when  no  experimental  rate  data  are  available. 
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